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We introduce three universality classes of chiral random matrix ensembles with a nonzero chemical 
potential and real, complex or quaternion real matrix elements. In the thermodynamic limit we find 
that the distribution of the eigenvalues in the complex plane does not depend on the Dyson index, 
and is given by the solution proposed by Stephanov. For a finite number of degrees of freedom, 
A'^, we find an accumulation of eigenvalues on the imaginary axis for real matrices, whereas for 
quaternion real matrices we find a depletion of eigenvalues in this domain. This eff'ect is of order 
1/^/N. In particular for the real case the resolvent shows a discontinuity of order l/\/iV. These 
results are in agreement with lattice QCD simulations with staggered fermions and recent instanton 
liquid simulations both for two colors and a nonzero chemical potential. 

PACS numbers: 12.38.Lg, ll.30.Qc, ll.15.Ha, ll.30.Rd 



Recently, nonhcrmitean random matrices with eigen- 
values scattered in the complex plane have received a 
great deal of attention in both condensed matter 
and QCD |^J^ . In condensed matter they have been used 
in problems ranging from neural networks [0 to the de- 
pinning transition of vortices in disordered superconduc- 
tors 1^,^. In QCD, they are relevant at nonzero chemical 
potential when the Dirac operator is nonhcrmitean. In 
this letter we introduce three universality classes of non- 
hcrmitean random matrices with the global symmetries 
of the QCD partition function. This allows us to interpret 
certain characteristic features of quenched Dirac spectra 
at nonzero chemical potential (i.e. with the fermion de- 
terminant ignored in generating the field configurations) 
as signatures of anti-unitary symmetries. Such signatures 
have been found in condensed matter applications as well 



In quenched simulations of QCD it appears that the 
critical chemical potential is proportional to the pion 
mass which vanishes in the chiral limit as {m is the 
quark mass) rather than to the baryon mass which re- 
mains nonzero for m ^ [^,^. This long-standing puz- 
zle has been resolved in a first successful application of 
random matrix theory to QCD at nonzero chemical po- 
tential The explanation is that the quenched 
limit is the limit Nf ^ of a partition function with 
the fermion determinant to the power Nf replaced by its 
absolute value. Such partition function develops a con- 
densate with Goldstone bosons consisting of a quark and 
a conjugate quark. Meanwhile, several other studies have 
confirmed this work 0. 

The above discussion was for Nc — 3. For Nc = 2 the 
situation is different . Then the Dirac operator is self- 
conjugate leading to a real fermion determinant. Because 
in this case also baryons can be Goldstone bosons, we 
expect a critical chemical potential of /ic ^ \/m. Math- 
ematically, the Dirac operator for Nc = 2 has an addi- 
tional anti-unitary symmetry. As is well-known in ran- 
dom matrix theory, the anti-unitary symmetries lead to 



three different universality classes [|i2|,|l3[ characterized 
by the so called Dyson index /?(/?= 1 for real, P — 2 
for complex, and /3 = 4 for quaternion real matrix ele- 
ments). In this letter we study a random matrix model 
of the QCD partition function at nonzero chemical po- 
tential for all three values of /3. 

The continuum Euclidean QCD partition function for 
Nf flavors with masses m / , and chemical potential fi can 
be written as an average over the Yang-Mills action, 



Nf 



Z{m, = ( det(7 ■ D + nif + fijo) 



(1) 



where 7 • is the Euclidean Dirac operator and 7^ are 
Euclidean Dirac matrices. In lattice QCD, the chemical 
potential is incorporated by including a factor in the 
forward time links and a factor e^^ in the backward time 
links. This implementation does not affect the symmetry 
relations discussed below. 

To obtain a random matrix theory corresponding to 
(|l|) we first write the Dirac operator in a chiral basis. 
Then the fi independent nonzero matrix elements are 
replaced by Gaussian distributed random variables cor- 
responding to the anti-unitary symmetries of the QCD 
partition function. In analogy with the random matrix 
partition function at zero ||l7| and nonzero temperature 
for Nf flavors, this partition function is defined as 
(see [fsl for a review) 



Z{m,n)= / DCP{C)Y[det{D{n) + 

f 



nif). 



(2) 



Here, C is an arbitrary N x N matrix with real, com- 
plex or quaternion real matrix elements. The integration 
measure DC is the Haar measure. The probability dis- 
tribution P{C) is given by 



P{C) = exp(-7V(/3/2)E2TrCC^) 



(3) 



1 



From now on we will work in units where E = 1. 

In the case oi Nc > 3 {P = 2) there are no anti-unitary 
symmetries and the Dirac operator is given by M 







iC + n 



iCT+n 



(4) 



Here, and in (6,8) below, /i is proportional to the identity. 

For Nc — 2, the Dirac operator in (|l|) is subject to an 
additional anti- unitary symmetry We stress that a 
nonzero chemical potential does not violate this symme- 
try. For continuum fermions (and naive lattice fermions 
pql) in the fundamental representation we have 



[CT2K,tjD{fi)]=0, 



(5) 



where C is the charge conjugation matrix, and i"fD{fi) 
is the continuum Dirac operator at ^ 0. Because 
(Cr2-ftr)^ = 1, it is possible to find a basis in which the 
Dirac operator is real. If we notice that is Hermitean 
we arrive at the following random matrix Dirac operator 







C + li 
H 



(6) 



where C is an arbitrary real matrix. Quenched instanton 
simulations at ^ are in this class . 

For staggered fermions the gamma matrices are absent 
and for Nc = 2 the anti-unitary symmetry is given by 



(7) 



where D^{ii) is the staggered Dirac operator at /i 7^ 0. 
Because {t2K)'^ = —1, we can organize the matrix ele- 
ments of the Dirac operator into real quaternions. For 
an (anti-)Hermitean Dirac operator the anti-unitary sym- 
metry results in a pairwise degeneracy of the eigenvalues. 
This degeneracy is broken at /i ^ 0. Then the eigenvalues 
of occur in complex conjugate pairs. In this symme- 
try class with /3 = 4 the random matrix model is given 
by 



C + fil 
-Ct + /il 



(8) 



where the matrix elements of C are real quaternions and 
1 is the unit quaternion. In addition to for Nc — 2, 
this class also contains the Dirac operator for gauge fields 
in the adjoint representation jis) for Nc> 2. 

In all three cases the spectral density at /x = is a 
semicircle which in the normalization defined by (H) (with 
S = 1) is given by 



Psc(A) - -^v/4-A2. 

ZTT 



(9) 



The above random matrix models apply to universality 
classes in which the eigenfunctions at /i = are extended. 



The eigenvalues corresponding to localized wave func- 
tions are statistically independent with spectral correla- 
tions given by the Poisson ensemble. A random matrix 
model in this class without a chiral structure was recently 
considered in [§J|]. At ^ 0, it was found that eigen- 
functions with complex eigenvalues are extended whereas 
eigenfunctions with real eigenvalues remain localized as 
for ^ = 0. The surprising result of this study was that 
a localization transition was observed in one dimension. 
In the present context the extended states can be in- 
terpreted as the emergence of a nonzero baryon number 
density when the eigenvalues scatter in the complex plane 
at /i 7^ 0. Another remarkable result of this study 0| 
was that in the thermodynamic limit a finite number of 
eigenvalues remained on the real axis for /z below a criti- 
cal value. Below we will show that the fraction of purely 
imaginary eigenvalues of the model (^ with [3 — 1 scales 
as ^ I/VN (Our convention and the conventions of 
for the eigenvalues differ by a factor i.) 

However, spectral correlations of the lattice QCD 
Dirac operator in four dimensions are given by the in- 
variant random matrix ensembles pO| ] even for relatively 
weak coupling. Moreover, explicit calculations of the in- 
verse participation ratio in an instanton liquid model 
shows delocalized Dirac eigenfunctions Therefore, 
it seems that the models in which the states are localized 
at fi — are inappropriate for QCD applications. 

In the remainder of this paper we study the spectrum 
of the random matrix model (|^) . The generating function 
is given by 



(det"(z--i:i)det"(z* -Dt)). 



(10) 



The resolvent is defined as G{z) = —dzlog Z/n and 
the spectral density in the complex plane is given by 
p(x,y) — dz-G[z)l'K. The quenched approximation is 
obtained as the limit n ^ at the end of the calculation. 
For integer n, this partition function can be analyzed in 
standard fashion [p2[P|. First we write the determinants 
as Grassmann integrals. Then we perform the average 
over the random matrices resulting in a 4-fermion in- 
teraction. If we bosonize this interaction and integrate 
out the fermions we arrive at a partition function that 
is amenable to a saddle point approximation. At /i 7^ 
this analysis was performed by Stephanov for Hermitean 
random matrices. The analysis for /3 = 1 and /3 = 4 
is somewhat more complicated ||2^. However, in both 
cases we succeeded to solve the saddle-point equations 
with the remarkable result that for the normalization de- 
fined by (^, the solution coincides with the one obtained 
for (3 = 2. For all values of /?, the resolvent can be eval- 
uated for arbitrary integer values of n. To leading order 
in 1/N we find the result is independent of n. Because of 
the absolute value of the determinant in ( p^ the gener- 
ating function is a smooth function of ri. This guarantees 
the existence of the replica limit. 
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FIG. 1. Scatter plot of the real (x), and the imaginary parts ( y) of the eigenvalues of the random matrix Dirac operator. 
The values of (3 and fj, are given in the labels of the figure. The full curve shows the analytical result for the boundary. 

ing from ten to one thousand. The full line represents 
the analytical result due to Efetov M 



Numerical results for the random matrix ensembles de- 
fined in are shown in Fig. 1. We show results for 
fi = 0.15 and ^ = 0.5. The dots represent the eigenvalues 
in the complex plane. The full line is the analytical re- 
sult 1^ for the boundary of the eigenvalues given by the 
algebraic curve (y^ -I- 4^^ ){fp — x^)"^ + = 4^^(^^ — a;^). 
For /3 = 1 and /3 = 4 we observe that the spectral den- 
sity deviates significantly from the saddle-point result. 
For /3 = 1 we find an accumulation of eigenvalues on the 
imaginary axis whereas for /3 = 4 we find a depletion of 
eigenvalues in this domain. This depletion can be under- 
stood as follows. For ^ = all eigenvalues are doubly 
degenerate. This degeneracy is broken at /.t 7^ which 
produces the observed repulsion of the eigenvalues. 

The number of purely imaginary eigenvalues appears 
to scale as VN, which explains that this effect is not 
visible in a leading order saddle point analysis. From a 
perturbative analysis of (|lO| ) one obtains a power series 
in 1/7V. Clearly, the V N dependence requires a truly 
nonperturbative analysis of Such a Vn scaling be- 
havior is typical for the regime of weak non-hermiticity 
first identified by Fyodorov et al. jsj. 

A similar cut below a cloud of eigenvalues was found 
in instanton liquid simulations for A'^c = 2 at ^ | |T^ ] 
and in a random matrix model of arbitrary real matrices 

. The depletion of the eigenvalues along the imaginary 
axis was observed earlier in lattice QCD simulations with 
staggered fermions p^]. 



In Fig. 2 we plot the fraction of imaginary eigenvalues 
versus /z-\/iV for ^ equal to 0.05, 0.1 and 0.2 and TV rang- 



PscWdX / dteM~N{2TTfipsc{X)?t% (11) 



ao 



where psc is defined in (^. This result was obtained 
for an ensemble of independent real symmetric matrices 
perturbed by an anti-symmetric matrix in the limit that 
the norm of the perturbing operator is of the order of the 
level spacing of the unperturbed matrix (his conventions 
differ from ours by a factor i). Apparently, the fraction 
of purely imaginary eigenvalues is not modified by the 
chiral structure of the ensemble (|^). Asymptotically, for 
^ 1, this fraction is given by ao l/p^/nN. 
It is well known that the replica trick fails in some cases 
p^ . For example, it fails for the unquenched partition 
function More typically, it fails in cases where the 
the saddle-point is given by a nontrivial manifold [ p4[ . 
In the limit /i — > and z real, the partition function 
has a higher degree of symmetry. Therefore, in the limit 
N ^ 00 with p'^N fixed, the solution of the saddle point 
equations is given by a nontrivial manifold, and there is 
no guarantee that the replica trick will work . In order 
to obtain truly nonperturbative results one has to rely 
on the supersymmetric method for random matrix the- 
ory |25[| . This method was extended to non-hermitean 
complex matrices in [^6|j|] and to arbitrary real matrices 
in The application of the supersymmetric method 
to the chiral ensembles will eventually provide us with 
an explanation of the scaling behavior of the number of 
purely imaginary eigenvalues ( |^7|] ). 
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FIG. 2. The fraction of real eigenvalues, ao versus 
for three different values of and matrices ranging from 
= 10 to iV = 1000. The full line shows Efetov's result 
[12]. 



Stephanov has shown analytically that the quenched 
approximation does not work for (5 = 2. Specifically, 
the unquenched partition function (||) for Nf ^ results 
in a nonzero chiral condensate below a critical value of 
fi. However in the quenched case, given by the n — > 
limit of (|l^), the chiral condensate is zero at any /i ^ 0. 
The situation for /3 = 1 and /? = 4 is different. Then 
in both cases the fermion determinant is real for real z. 
(For /3 = 4 this follows by using the identity q* = a2qcr2 
for a quaternion real element q.) Therefore, for real 2, 
the generating function for n replicas is identical to 
the partition function (||) for 2n flavors with mass z. As 
mentioned above, the saddle point result for the resolvent 
defined by the generating function (^0|) does not depend 
on n and gives the quenched result for n 0. We thus 
conclude that quenching works for an even number of 
flavors. Consequently, chiral symmetry will be restored 
for arbitrarily small nonzero /x, whereas a condensate of 
a quark and a conjugate quark develops. Indeed, this 
phenomenon has been observed in the strong coupling 
limit of lattice QCD with two colors jl^,^ . 

In conclusion, at ^ we have found that the deple- 
tion of Dirac eigenvalues on the imaginary axis observed 
in lattice QCD simulations and the accumulation of Dirac 
eigenvalues found in instanton liquid simulations (both 
for two colors and quenched) is a generic feature related 
to the anti-unitary symmetries of the Dirac operator. If 
it turns out that the Dirac eigenfunctions are localized 
at (U = (see ||2^ for lattice results in this direction) the 
fraction of purely imaginary eigenvalues at 7^ might 
remain flnite in the thermodynamic limit. For two col- 
ors, this might lead to unexpected chiral properties in 
the continuum limit at /i 7^ 0. Clearly, this possibility 
deserves further attention. 
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